Reconciling motion and stationarity in single-photon polarization states 
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The question of the representation of quantum stationary partially polarized waves as random 
superpositions of different polarization ellipses is addressed. To this end, the Bohmian formulation of 
quantum mechanics is considered and extended to quantum optical polarization. As is shown, this 
approach properly combines definite time evolving trajectories with rigorous stationary quantum 
distributions via the topology displayed by the associated phase field. 
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I. INTRODUCTION 

In most courses, textbooks, and specialized treatises 
it is common to introduce the polarization of harmonic 
waves as an ellipse described by the electric field in a real 
configuration representation [l|. Within this approach, 
partial polarization thus arises by the rapid and random 
succession of more or less diff'erent polarization states. In 
more advanced approaches, though, the Stokes parame- 
ters and Poincare sphere are introduced, which allow us 
to assess the degree of polarization and make use of the 
powerful tools provided by the SU(2) group. Moreover, 
the Stokes parameters also enable a simple connection to 
the quantum regime in terms of their quantum counter- 
parts, namely the Stokes operators [2]. 

Once the Stokes parameters are introduced, usually the 
polarization ellipses are left behind, so that the statistics 
of the electric field in the polarization plane are aban- 
doned. This is particularly remarkable within the quan- 
tum domain, where uncertainty relations imply that no 
field state can describe a perfect ellipse, just in the same 
way that no particle can follow a definite trajectory 043 ■ 
This thus suggests some kind of quantum mismatching 
between the intuitive ellipse picture and the more power- 
ful Stokes-Poincare formalism. More specifically, we find 
quantum states with maximum degree of polarization, 
but with their electric-field distribution far from resem- 
bling any ellipse, as seen, for example, in Fig. [5] below 
(for other examples, see also Fig. 1 in Ref. y or Fig. 4 in 
Ref.|j). Furthermore, the mismatch seems aggravated for 
quantum stationary states (i.e., states with definite total 
photon number). In such a case, the electric-field distri- 
bution is constant even for states with maximum degree 
of polarization, which should correspond to an electric 
field rapidly describing an elliptical trajectory. 

In this work we investigate whether there is still a pos- 
sibility to keep the most intuitive approach, where par- 
tially polarized waves are devised as random superposi- 
tions of different polarization ellipses. To this end, the 
Bohmian formulation of quantum mechanics loj or, more 
properly speaking, its extension to optics [7|, [Sj , seem to 
be very convenient. This approach properly combines the 
two ideas that we want to combine: definite trajectories 



and rigorous quantum distributions. More specifically, 
here we are going to study the Bohmian trajectories (or 
optical paths) for a two-dimensional (2D) isotropic har- 
monic oscillator, which properly represents in the quan- 
tum domain a transversal two-mode harmonic wave. For 
definiteness, and for the sake of simple illustration, we 
consider single-photon pure states. Accordingly, we have 
organized this work as follows. In Sec.|lT]we introduce the 
basic elements involved in the representation of single- 
photon polarization states. In Sec. IIIII the associated 
Bohmian dynamics is analyzed and discussed. Finally, 
in Sec. IIVI the main conclusions drawn from this work 
are summarized. 



II. POLARIZATION OF ONE-PHOTON STATES 

A. Polarization ellipse 

Consider a harmonic light wave consisting of two 
modes of the same frequency, w, and with their corre- 
sponding electric fields vibrating along orthogonal direc- 
tions. These two modes are represented by their complex 
amplitude operators ai and 0,2- This system is equivalent 
to a particle in a 2D isotropic harmonic potential. The 
equivalence becomes clearer through the quadrature op- 
erators representing the real and imaginary parts of the 
electric field, a^ 



-] oc Xj 



iYj, with 



A, 



1 

V2 



{a]+a,y 



K,- 



V2 



(«] - a.) , (1) 



and commutator [Xj, Yk] = iSjk, with j,k = 1,2. The Xj 
and Yj operators are thus formally equivalent, respec- 
tively, to the (dimensionless) position and momentum 
of a 2D massive particle. For simplicity, and to exploit 
as much as possible this equivalence, we assume h — 1, 
uj = 1, and m = 1, with the latter being the mass of the 
equivalent effective particle. The corresponding effective 
Hamiltonian then reads as T-L = ajoi -f 0202. 

The quadrature operators, Xj (with j — 1,2), also 
allow us to introduce a wavefunction for the 2D real 
electric field by projection of the field state jf/j) on the 



unnormalized joint eigenstate of the quadrature opera- 
tors Xj\xi,X2) = Xj\xi,X2)- In particular, for photon- 
number eigenstates, with a-aj\ni,n2) — nj\ni,n2), we 
have 

(a;i,X2|ni,n2) =A/'i2i?„i(xi)iJ„,(a;2)e-(^'+"')/2, (2) 

where A/12 = (2"i+"27rni!n2!)~^/^ is the norm and Hn 
are the corresponding Hermite polynomials. Accordingly, 
the quantum analog of the polarization ellipse is the dis- 
tribution for the 2D real electric field [^-Q, 



p{xi,x2,t)^\{xi,x2mm^, 

with IV'(i)) = exp(-im)|V'(0)). 

B. Stokes parameters 



Regarding the Stokes picture of polarization, it can be 
well started from the Stokes operators [4I, 



So - 

Sy — 

s. = 



a\ai 
a\ai 



a\a2, 
a\a2, 
a2a\, 
i{d2a[ — d^ixi) 



which satisfy the commutation relations 
with cyclic permutations, and 



[Sx^ Sy\ — 



S2 = 5'o(5'o + 2I), 



[S,^o]=6, 



(4) 



(5) 



(6) 



with S ~ {Sx, Sy, Sz)- The classical Stokes parameters 
are the mean values of the Stokes operators, Sj = (Sj). 
Because of the non-vanishing commutator, as described 
by Eq. (O, no state can have definite values of all Stokes 
operators simultaneously (the only exception being the 
two-mode vacuum, where they vanish trivially). This is 
conveniently expressed by the uncertainty relation 



(ASf = (ASx)^ + (ASyf + [ASzf 
= (^0(^0 + 2i))-(S)2> 2(^0), 



(7) 



which holds after Eq. (jH]) 9] taking into account that 

KS)| < (So). 

Consider now the standard (classical) definition of the 
degree of polarization is 7^ = |s|/so- Some other more 
complete definitions, though, have also been proposed, 
particularly within the quantum domain [J, |5|, |l2| . Based 
on the fact that \s\ < sq, the Stokes parameters readily 
provide a representation for polarization states in a unit 
sphere, namely the Poincare sphere. This is done through 
the parametrization 



Sx = socos(2x)cos(293), 
Sy = socos(2x)sin(2<^), 



sosin(2x). 



(8) 





(3) FIG. 1: Polarization ellipse (left) and Poincare sphere (right). 



which is sketched in Fig. [T] In the transformation re- 
lations ([5|), tan(2(^) — Sy/sx, with ip being the angle 
between the major axis of the polarization ellipse and 
the xi axis. Regarding x, this angle determines the ratio 
between the minor and major axes of the polarization el- 
lipse (6 and a, respectively), with 7r/4 > x ^ — 7r/4 and 
tanx = ±6/a, where the sign is given by the handedness 
(right-handed if Sz > and left-handed if Sz < 0). The 
relations (|8]) thus provide us with a one-to-one correspon- 
dence between the Stokes parameters (or points on the 
Poincare sphere) and some average or mean polarization 
ellipse. 



C. One-photon states 

Let us now consider, more specifically, the most gen- 
eral pure one-photon state. In the photon-number basis 
\ni,n2) it reads as 



|^)=Ci|l,0)+C2|0,l), 

with Stokes parameters 



50 = |C1|2-H|C2|2 = 1, 
Sx = \C1? - C2 P, 



Sy — c*C2 -I- C1C2 = 2Re(ciC2), 
Sz = i(cjc2 — C1C2) = 2Im(ciC2). 



(9) 



(10) 



These are all stationary states, i.e., 'H\4') = IV'): so that 
\il}{t)) = exp(-ii)|V'(0)). They are also SU(2) coherent 
states 13], which are usually considered as classical-like 
regarding polarization [14| . as well as being minimum 
uncertainty states of the uncertainty relation ([7|) . More- 
over, since \s\ = sq for these states, we have 7-" = 1, thus 
displaying the maximum degree of polarization accord- 
ing to the classic definition seen above. In spite of this, 
they cannot be considered as having perfect polarization, 
since (AS)^ — 2(6*0) 7^ 0. This is properly reflected by 
more complete assessments of the degree of polarization 

ai[ii. 

The quadrature wavefunction associated with the pure 
one-photon state ^ is 



{xi,X2\lp(,t)) 



{cixi + C2X2) e 



-{xl+xl)/2-it 



(11) 




D ~^^^^^^^^^^^^^^^^^^^^^^^^^ -^^^^^^^^^^^^^^^^^^^^^^^^^ ~^^^^^^^^^^^^^^^^^^^^^^^^^ 



the equator of the sphere this is for hnearly polarized 
Ught). In this case, P{xi,X2) essentiaUy consists of two 
separate lobes. For 2if — these lobes are parallel to 
the xi-axis (see lowest left panel), but their inclination 
changes as 2ip increases, thus becoming parallel to the 2:2- 
axis when 2ip — t: (for the intermediate case, 2ip — 7r/2, 
they are 7r/4-inclined, as seen in the lowest right panel). 
Apart from deviating from the ellipsoidal picture, it is 
also worth stressing that information about the handed- 
ness will also get lost. This is because P{xi,x2) does not 
depend on Sz- 



III. BOHMIAN DYNAMICS 

A. Equations of motion 

Let us now consider the Bohmian approach with the 
purpose of grasping some physical insight from the dy- 
namical state associated with the wavefunction (|lip . 
Within this context, the corresponding Bohmian trajec- 
tories, r{t) = [xi{t),X2{t)], are obtained after integration 
of the guidance equation 



-3-2-101 2 3 



f = V5 = 



-ie^'^Ve*^, 



(13) 



FIG. 2: (Color online) Contour-plots corresponding to the 
P{xi,X2) distribution of the 2D real electric field associated 
with the pure one- photon state ([9]). The values chosen for 
the angles 2x and 2ip cover the first quadrant of the Poincare 
sphere (see Fig.[l] right). Minimum to maximum is indicated 
with the transition from blue to red, starting from zero in all 
cases. 



where S is the (real-valued) phase of the wavefunction 
(HH), i.e.. 



js 



(CiXi + C2X2)e 



^\cxYx\ -H \C2\^X\ + (CIC^ + C|C2)X1X2 



(14) 



In terms of the Stokes paranreters, this expression can 
also be recast as 



and the corresponding electric-field distribution is 



P{xy,X2) 

-[(1 + S.)r2 
TT 



[(1 -f s^)x\ + (1 - s,)x2 + 2syX^X2\ e-(-?+-'), (12) 



which is effectively independent of time, i.e., it is station- 
ary, as expected. In Fig. [2] we show a series of contour- 
plots of P{xi,X2) for angles 2% and 2(p that cover the 
first quadrant of the Poincare sphere; the behavior in the 
remaining quadrants is analogous and can be inferred 
by taking into account the corresponding synrmetry re- 
lations. Also, for convenience regarding the computation 
of the Bohmian trajectories (see discussion below), we 
have considered the minimum 2% to be tt/IO instead of 
0. As can be seen, only for 2% — ±7r/2 and any (p the 
picture of an ellipse remains valid (this is for circularly 
polarized light). As the value of these angles is changed, 
this picture breaks down despite having the maximum 
degree of polarization in any of them {V = 1), with the 
most dramatic breaking-off being for 2% = 0, i.e., along 



VT 



JS 



Xi 



1 



X2 



-zi+i(5 



^yxl + x1+ Sx{x\ - X^) + 2SyXlX2 



(15) 



where i5 is a global relative phase associated with the 
coefficient ci . This phase factor is physically meaningless 
regarding the topology displayed by the phase field, S, as 
well as the trajectory dynamics described by the equation 
of motion below, as also happens with t. Substituting 
this expression into Eq. (|13l) . we find 

{X2,-Xl). (16) 



CI+ xl+ Sj; {xl — xl) + 2SyXlX2 



In general, Eq. ([TS]) has to be numerically integrated 
in order to obtain the corresponding trajectories, as seen 
in Sec. nil Cj below. However, as shown in next Section, it 
is also possible to draw a series of interesting conclusions 
directly from the form of Eq. (|16l) . without the need of 
integrating it. 



B. Properties of the trajectories 

1. Circular trajectories 

In spite of the complex dependence on Xi and X2 dis- 
played by the prefactor of the equation of motion (1161) . 
the trajectories for any single-photon state with Sy ^ Q 
are always circular. This comes from the fact that 



VS* 



0, 



(17) 



which implies that r^ (and therefore r) is constant. 

Now, for Sz — 0, when the photons are linearly po- 
larized, we have r = and therefore the whole vector 
r will be constant. That is, each separate r-component 
is constant with time, since Eq. (1171) only ensures the 
time-independence of r — |r|, but not its components 
(see below). This case can be then regarded as the limit 
of a circle described at a vanishing speed. 



2. Nonuniform angular frequency 

The angular frequency of these circular trajectories is 
not uniform. This is readily seen if we use polar coordi- 
nates, 



Xi 



r cos< 



X2 



(18) 



and then express ([TC)l in terms of these coordinates, 
which renders the following equation of motion for the 
polar component 



Sz/r'^ 



\ + Sx COS 2(j) + Sy sin 20 



sin(2x)/r^ 



1 + cos(2x) cos(20 - 2ip) ■ 



(19) 



We recall that the polarization ellipse is described at con- 
stant angular speed uj (in our case uj — 1). However, as 
seen through ^9]) . this is not the case for the Bohmian 
trajectories, where cf) strongly depends on the polariza- 
tion state (see below). Thus, the only case with uniform 
angular velocity holds for x = i'''/4, i- e., for circularly 
polarized light. 



3. Handedness 

The handedness of the Bohmian trajectories is the 
same than for the corresponding polarization ellipses. Af- 
ter the change of coordinates p^ . the Bohmian motion 
will be right-handed if < 0, while the polarization el- 
lipse is right-handed when s^ > 0. 



4- Extreme instantaneous frequencies 



and ellipse parameters 

The maximum instantaneous angular frequency of the 
Bohmian trajectories, denoted by il = \(f)\, is 



^^max — 



1 



r2|tanxr 



(20) 



which holds when cos(20 — 2(p) — —1, i.e., when </> = 
if ± 7r/2, modulo n. In other words, this condition is 
equivalent to say that the maximum instantaneous an- 
gular frequency of the Bohmian trajectory coincides with 
the direction of the minor axis of the polarization ellipse 
(see Fig. [U. 

On the contrary, the minimum instantaneous angular 
frequency is 



n. 



|tanx| 



(21) 



which holds when cos(2(/) — 2ip) = 1, i.e., when cj) = ip, 
modulo TT. That is, the minimum instantaneous angular 
frequency coincides with the direction of the major axis 
of the polarization ellipse. Note, therefore, that the ratio 
between minimum and maximum instantaneous angular 
frequencies is 



^n 



n^ 



tan X = ^, 



(22) 



which coincides with the ratio of the minimum to max- 
imum axes of the polarization ellipse, b and a, respec- 
tively. 



C. Numerical simulations 

A priori, the properties discussed above may seem 
counterintuitive and puzzling. For example, how is it 
possible that the trajectories describe circles if their an- 
gular frequency is variable with time, or if they tend to 
distribute in regions where P(xi, 2:2) is maximum, avoid- 
ing the nodal ones? Actually, it seems that there is a 
mismatching between the symmetry of the wavefunction 
and the features characterizing the trajectories. In order 
to elucidate these questions as well as to reconcile two an- 
tagonists, namely motion and stationarity in the case of 
single-photon states accounted for by wavefunctions like 
([9]) , we have numerically integrated the equation of mo- 
tion (|16p for some of the conditions considered in Fig. [21 
Each condition represents a different single-photon state 
vector and, therefore, a particular dynamics. The corre- 
sponding trajectories as well as some other additional 
representations of interest are displayed in Figs. [31 [H 
and [5l For clarity, the initial condition, xi (0) — 1 and 
2^2(0) = 0, was chose to be the same in the three cases. 
Accordingly, all the trajectories will be circles of unit 
radius and, in principle, indistinguishable, in agreement 
with (|17p . Disambiguation follows, though, when one an- 
alyzes the time-dependence of each component separately 



or their respective velocities, as can be seen in the addi- 
tional representations (see also discussions below). For 
completeness, and also to verify our conclusions in other 
cases, we also considered the same cases, but with other 
radii. The corresponding trajectories will not be plotted 
here, although we would like to stress the fact that we 
could corroborate that, effectively, the smaller the radius, 
the faster the motion around the circle, in agreement with 
the r~^-factor in Eq. (fT9|) . Moreover, it has also been no- 
ticed that the time needed for trajectories with the same 
radius to complete the circle is different, increasing as 2% 
approaches zero. This can be seen below by examining 
the dependence on time of the velocity components of the 
trajectories, i.e., Vi = i^, with i = 1,2 [see red dashed 
line in panels (c) and (d) of Figs. [31 IH and[S]. Thus, 
as 2x approaches zero, one observes a certain bistable 
behavior: vi and V2 undergo important variations in ex- 
tremely short periods of time, while they remain finite (or 
even vanishing) for relatively long times, thus conferring 
certain stationarity or stability to the corresponding pho- 
ton state. This fact is actually in correspondence with 
the last of the properties studied in Sec. IIIIBl 

Going now to each particular case, in Fig. [3] we have 
displayed the results for circularly polarized light 2x — 
tt/2 and 2(p — on thcPoincare sphere (see Fig[T|). Notice 
that the evolution is ruled by a harmonic motion. This 
can be seen through the time-dependence of the compo- 
nents xi and X2, and their corresponding velocities vi 
and V2 [see Figs.jljc) and (d)], or through the respective 
phase-space orbits [see Figs. [3l^e) and (f)]. In this case, 
Sz = So = 1 and s^ ^ Sy = 0, which means that the pho- 
ton has equal probability to be in the |l,0)-state than 
in the |0, l)-state. Therefore, this uncertainty leads the 
trajectory to visit all points on a circle of a given radius 
at the same (angular) velocity, this eventually manifest- 
ing as a harmonic-like motion or oscillation between 1 1 , 0) 
and |0, 1) (or, in other terms, between the horizontal and 
vertical polarization states). 

As we move towards the equator of the Poincare 
sphere, i.e., vanishing Sz for linearly polarized light (from 
top to bottom in Fig. [5]), the torus- like distribution starts 
developing two lobes, which in the limit Sz — become 
separate. Here we find an apparently paradoxical behav- 
ior: P{xi,X2) consists of two separate lobes, but the tra- 
jectories are still circles, as indicated by ([T7)) . To analyze 
this situation, we have proceeded as before, but consid- 
ering the single-photon state defined by 2% = tt/IO and 
2ip = (see Fig.|3]). In order to reconcile both behaviors, 
principally because no trajectory should be expected at 
a nodal region (these are regions of lowest probability), 
it is important to observe the time-dependence of xi and 
X2, displayed in Figs. Hfc) and (d), respectively. In the 
case of xi, we note a bistable behavior. That is, xi is 
essentially -1-1 or —1 for relatively long periods of time 
(which are expected to increase as 2^ — >■ 0), while X2 only 
oscillates linearly, up-and-down, between these values to 
compensate the transition from Xi to X2- We observe 
that such transitions are extremely fast. Actually, in 
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FIG. 3: (Color online) Different representations of the 
Bohmian trajectories associated with a single-photon polar- 
ization state characterized by 2x = 7r/2 and 2(fi — on the 
Poincare sphere. In panels (c) and (d) the black solid lines de- 
notes the spatial component {xi or X2), while the red dashed 
line refers to the respective velocity component {vi or 112). In 
the corresponding panels, the velocity components have been 
obtained evaluating (|16|) along the trajectory. 



the case of xi, the velocity vi is negligible except at the 
transition times, while V2 fluctuates between two almost 
stationary values (-1-0.16 and —0.16, approximately). As 
a consequence, the system will remain apparently either 
in one lobe or the other of P{xi,X2), as resembled by the 
orbit of the velocity representation displayed in Fig. SJ^b) 
(if we consider 2ip — ?> it, we would find the same, but 
exchanging xi by X2 and vice versa). 

The situation is similar if we move around the equa- 
torial plane regarding the presence of the two lobes in 
P{xi, X2), although this motion implies a rotation of the 
axis along which they are distributed. As mentioned 
in Sec. Ill CI the lobes are aligned with the xi-axis for 
2ip = and start get inclined anticlockwise as this an- 
gle increases, becoming aligned with the a::2-axis when 
2(p = TT. This change does not only affect the distribu- 
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FIG. 4: (Color online) As in Fig. [1 but for 2x = tt/IO and 
2(p = 0. 



FIG. 5: (Color online) As in Fig. El but for 2x = tt/IO and 

2ip = 7r/2. 



tion P(xi,X2), but also the phase of the single-photon 
state and therefore its associated trajectory dynamics. 
To illustrate this case, we are going to consider the state 
described by 2% — tt/IO and 2(p — tt/2 (see Fig. [5]), which 
is linked to the distribution displayed in the lowest right 
panel of Fig. [2] Although the trajectory is still circu- 
lar, the orbit in the velocity space describes two lobes 
distributed along a 7r/4-axis [see Fig. EJb)], which in 
the limit 2% — ^ will end up in a line with a 7r/4- 
inclination. Note that these lobes constitute a distor- 
tion of those displayed by the velocities in the previous 
case [compared with Fig. |4l[b)]. In Figs.[5l[c) and (d) we 
observe that this behavior just means that the photon 
oscillates between one lobe of the distribution P{xi,X2) 
and the other in a rather inhomogeneous way: there is 
a slow, gradual approach to the transition followed by 
a sudden switch, as indicated by the corresponding ve- 
locities. In the corresponding phase-spaces [see Figs.[5][e) 
and (f)], we observe precisely that the velocity undergoes 
its maximum value right after the photon has reached one 
of the lobes [of P{xi,X2)] and just before it jumps into 



the other. Nonetheless, although the angular velocity is 
rather nonuniform, the trajectory is still circular. 

Finally, we would like to illustrate the source of motion 
in the kind of stationary states considered here. Con- 
trary to non-stationary wavefunctions, where the time- 
dependence of the phase leads naturally to time-evolving 
trajectories [H, [Taj , here motion has a topological ori- 
gin: the space-dependent (configuration space) gradient 
of the global phase associated with the stationary state. 
In Fig.[S]we show the distribution P{xi, X2) (upper row), 
the phase S (second row), and the components of the lo- 
cal velocity field, vi and V2 (in the third and fourth rows, 
respectively), of the three cases represented in Figs. [3l 
m and [S] (left, central, and right columns, respectively). 
From a quick inspection, the role of the phase field topol- 
ogy on the trajectory dynamics becomes readily appar- 
ent. This field is decreases monotonically counterclock- 
wise, thus explaining the evolution of the trajectories 
(from their starting, as xi decreases, X2 increases). Now, 
while this decrease is gradual for 2% = 7r/2, it becomes 
step-like as 2% approaches zero. In the latter case, the 
step lies precisely along the symmetry axis passing be- 
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FIG. 6: (Color online) Distribution P{xi,X2) (upper row), 
phase field S (second row), and components of the local ve- 
locity field, vi and V2 (third and fourth rows, respectively), 
for the polarization states considered in Figs.[3l|4l and[5](left, 
central, and right columns, respectively). Minimum to max- 
imum is indicated with the transition from blue to red: for 
P{xi,X2), starting from zero in all cases; for S, between —tt 
and +tt; and for vi and V2, between —2 and +2. In the last 
two rows the shaded areas indicate a truncation to the maxi- 
mum and minimum values chosen. 



tween the two lobes. In this sense, the motion is rel- 
atively slow along each plateau of the phase field and 
very fast when the trajectory conies down the step (i.e., 
passes from the region covered by one of the lobes to 
the region covered by the other). This has a straight- 
forward counterpart in the corresponding velocity field 
components. For 2% = these components are symmet- 
ric with respect to a 7r/2-rotation, which is a signature 
of the harmonic motion (as one of the components starts 
to decrease, the other increases, and vice versa). Now, as 
2x approaches zero we find two types of behaviors. For 
2x = and 2(p — 0, the component along the step of the 
phase field becomes more prominent than the component 
that is perpendicular to this step, which eventually be- 
comes meaningless. This is in correspondence with the 
fact that it is more likely to find the trajectory either 
in the right or the left lobe of P{xi,X2) (small vi and 



V2), but not in between (large vi). A similar behavior, 
but exchanging the axes would be found for 2^ = with 
2(y3 = vr. Now, as 2ip moves towards 7r/2, both compo- 
nents tend to align, thus becoming equal for 2ip — 7r/2, 
which means that the trajectories will avoid to stay in a 
neighborhood of the diagonal X2 = —xi. 



IV. CONCLUSIONS 

In this work we have shown that the topology displayed 
by the Bohmian trajectories analyzed have no relation 
whatsoever with the form of the polarization ellipse of 
the corresponding field state, since they are always cir- 
cles. However, the angular speed at which these circles 
are described contains complete information about the 
mean polarization state, i.e., the ellipse traditionally as- 
sociated with the Stokes parameters. This information 
is paradoxically absent from the probability distribution 
corresponding to the 2D real electric field [see Eq. ([T2|) ]. 
where one would expect to see the quantum counterpart 
of the classical mean ellipse. Moreover, seemingly, the 
Bohmian trajectories provide no hint about the uncer- 
tainty of the polarization state, since the family of tra- 
jectories associated with a given field state are essentially 
identical. 

Concerning the dynamics associated with single- 
photon polarization states, we have found a relatively rich 
variety of behaviors depending the point on the Poincare 
sphere defining the state vector. As seen, although all 
the trajectories are circles, their dynamical evolution is 
strongly dependent on their position with respect to the 
Poincare sphere, ranging thus from harmonic behaviors 
(for trajectories on the poles) to bistable oscillating ones 
(for trajectories on the equator). In all cases, though, 
this dynamics has a clear topological origin associated 
with the particular shape of the phase field, which goes 
from a gradual decrease around (a;i,X2) = (0,0) (at the 
poles) to a step- like one (on the equator). Accordingly, 
trajectories evolve cotmterclockwise either harmonically 
or jumping-like, respectively. 

The fact that polarization stationary states may have 
an associated internal or intrinsic dynamics is puzzling 
as well as challenging, for it goes beyond our physical 
intuition even if it is proven mathematically, as done 
here. Even though distributions like P(xi,X2) are in- 
variant objects with time, there is a dynamics induced by 
the phase field topology. This posses the interesting and 
stimulating question of whether such an intrinsic motion 
could be experimentally measured. Note that recently 
the group of Steinberg at the University of Toronto was 
able to experimentally infer [131 the Bohmian trajecto- 
ries or averaged photon paths [7| in the two-slit experi- 
ment by means of weak measurements [18|,[l3- The ideas 
posed here clear point towards a further extension of this 
work in this direction with the ultimate goal to perform 
such experimental measurements. This would allow us 
to better understand the nature of stationary states as 



situations where the system is statisticaUy described by 
a time-independent (steady) distribution, but that dis- 
play particular inner dynamics. Obviously, this has a 
potential interest in quantum information and quantum 
computation, where the basic ingredient, the qubit, is 
precisely described by a state like ([9]) . Furthermore, also 
notice the implication in quantum mechanics, where the 
role equivalent to polarization is the particle spin. This 
vectorial quantity is traditionally assigned to an inter- 
nal rotation, because it fulfils the same properties of the 
rotations group. With the picture on stationary states 
provided above, based on a dynamical ground, this con- 
nection becomes even clearer. 
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